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Genetic factors impact liver injury susceptibility 
and disease progression. Prominent histological 
features of some chronic human liver diseases 
are hepatocyte ballooning and Mallory-Denk bodies. In 
mice, these features are induced by 3,5-diethoxycarbonyl- 
1 ,4-dihydrocollidine (DDC) in a strain-dependent man- 
ner, with the C57BL and C3H strains showing high and 
low susceptibility, respectively. To identify modifiers of 
DDC-induced liver injury, we compared C57BL and C3H 
mice using proteomic, biochemical, and cell biologi- 
cal tools. DDC elevated reactive oxygen species (ROS) 
and oxidative stress enzymes preferentially in C57BL 
livers and isolated hepatocytes. C57BL livers and he- 
patocytes also manifested significant down-regulation, 



aggregation, and nuclear translocation of glyceralde- 
hyde 3-phosphate dehydrogenase (GAPDH). GAPDH 
knockdown depleted bioenergetic and antioxidant en- 
zymes and elevated hepatocyte ROS, whereas GAPDH 
overexpression decreased hepatocyte ROS. On the other 
hand, C3H livers had higher expression and activity 
of the energy-generating nucleoside-diphosphate kinase 
(NDPK), and knockdown of hepatocyte NDPK aug- 
mented DDC-induced ROS formation. Consistent with 
these findings, cirrhotic, but not normal, human livers 
contained GAPDH aggregates and NDPK complexes. 
We propose that GAPDH and NDPK are genetic modi- 
fiers of murine DDC-induced liver injury and potentially 
human liver disease. 



Introduction 

Race and ethnicity have been shown to be important factors in 
the development and progression of liver diseases for reasons 
that go beyond socioeconomic status (Nguyen and Thuluvath, 
2008). Understanding the biological basis for susceptibility to 
liver diseases may provide insights into novel therapeutic strat- 
egies for their prevention and treatment. Alcoholic liver disease 
(ALD) and nonalcoholic steatohepatitis (NASH), which repre- 
sent a significant societal burden, also have a strong genetic 
component (Anstee et al., 2011). Hepatocyte ballooning and 
protein aggregation in the form of Mallory-Denk bodies 
(MDBs) are common characteristics of the liver injury pheno- 
type and correlate with worse outcomes in patients with ALD 
and NASH (Cortez-Pinto et al., 2003; Gramlich et al., 2004). 
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However, MDB accumulation is not observed in all patients 
with the same liver disease (Gerber et al., 1973; Fleming and 
McGee, 1984), and MDBs may be more prevalent in Hispanics 
in comparison with other groups (Mohanty et al., 2009; Rakoski 
etal.,2011). 

One approach toward understanding the biological basis 
for these genetic differences in patient populations is to delin- 
eate the mechanisms of hepatocyte injury in animal models 
with varying genetic susceptibility. To that end, we previously 
found that experimental induction of hepatocyte ballooning and 
MDBs in mouse livers by using an established model of chronic 
administration of 3,5-diethoxycarbonyl-l,4-dihydrocollidine 
(DDC; Yokoo et al., 1982; Zatloukal et al., 2007) manifests in a 
highly strain-dependent manner (Hanada et al., 2008). Compar- 
ison of five inbred mouse strains (Beck et al., 2000) with 
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Figure 1 . Proteomic comparison of C3H and C57BL mouse livers before and after DDC treatment by 2D DIGE analysis and classification of differentially 
expressed proteins into function-based groups. (A) Liver homogenates were prepared from untreated mice (Control) or mice fed a DDC-containing diet for 
3 mo. The liver proteins from the C3H and C57BL mice were labeled with the fluorescent dyes Cy3 (green) and Cy5 (red), respectively. The circled spots 
indicate differentially expressed proteins that were selected for identification by MS. (B) From an initial analysis of 80 selected gel spots, proteins that were 
identified as having significant expression differences (ratio greater than two) between the two mouse strains under basal conditions or after DDC treatment 
are shown. Note that PRDX6 is represented twice, reflecting the acidic and basic forms of the protein. (C) The proteins shown in B were classified into three 
groups based on their known functional properties. 



distinct genealogies (FVB/N, C3H/He, Balb/cAnN, C57BL/6, 
and 129Xl/Sv) for their susceptibility to DDC-induced hepato- 
cyte ballooning and MDBs identified the C57BL/6 and C3H/He 
strains (hereafter referred to as C57BL and C3H) as being most 
and least susceptible, respectively (Hanada et al., 2008). The 
molecular constituents of MDBs, which mainly include cyto- 
skeletal and stress proteins, indicate that they share common char- 
acteristics with nonepithelial cell inclusions of inherited and 
age-related brain and muscular diseases (Muchowski and Wacker, 
2005; Askanas et al., 2009; Goldfarb and Dalakas, 2009; Liem 
and Messing, 2009). Therefore, elucidating the factors that 
likely underlie MDB formation in hepatocytes may also provide 
insights into neuro- and myodegenerative disorders involving 
protein aggregation. 

In the current study, we performed 2D differential in-gel 
electrophoresis (DIGE) and mass spectrometric proteomic com- 
parison of healthy and DDC-injured livers of MDB-susceptible 
(C57BL) and MDB-resistant (C3H) mice. We then evaluated 
the differentially expressed proteins and their corresponding 
signaling pathways in isolated hepatocytes and in normal and 
diseased mouse and human livers. Our findings indicate that ele- 
vated oxidative stress impairs energy metabolism by promoting 



glyceraldehyde 3-phosphate dehydrogenase (GAPDH) down- 
regulation and redistribution from a soluble cytoplasmic form 
into insoluble aggregates residing in the cytoplasm and nucleus. 
Furthermore, we show that GAPDH is an upstream regulator of 
several antioxidant and energy metabolism enzymes, including 
nucleoside-diphosphate kinase (NDPK). Therefore, inadequate 
supply of energy equivalents likely diminishes the cell's capac- 
ity to properly fold and degrade proteins, likely leading to hepa- 
tocyte injury and protein aggregation in the form of MDBs. 

Results 

Proteomic comparison of C3H and 
C57BL mouse livers before and after 
DDC treatment 

A 2D DIGE analysis was conducted to compare the liver pro- 
teome of C3H (MDB resistant) and C57BL (MDB susceptible) 
mice. This approach yielded a large number of differentially 
expressed proteins in untreated and DDC-treated mice from the 
two strains (Fig. 1 A). The circled green- and red-colored spots 
(Fig. 1 A) were subsequently subjected to tandem mass spec- 
trometry (MS/MS) for protein identification and determination 
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Figure 2. Strain-specific and DDC-induced changes in the expression of regulators of protein processing, energy metabolism, and oxidative stress. 

(A-C, top panels] Equal amounts of total liver protein from three untreated and four DDC-fed C3H or C57BL mice were resolved on SDS-PAGE gels and 
immunoblotted with antibodies to the indicated proteins. Blots separated by solid lines indicate nonconsecutive lanes of the same equally exposed mem- 
brane. The expression of some proteins (PDIA4, FAH, and ME1) was relatively unaffected across groups, and they serve as loading controls, (bottom 
panels) Relative band intensity values were plotted for each protein in the three different groups. ME1 , malic enzyme 1 ; mono, monomer; oligo, oligomer. 
*, P < 0.05; **, P < 0.01 ; ***, P < 0.001 using a two-way analysis of variance. Each tested group included three to four mice. Results are represented 
as the mean and the SD. 



of the relative expression ratios between the two strains. Of 80 
targets whose identity was characterized by MS/MS, those that had 
the highest expression difference either before of after DDC 
treatment were selected for further characterization (Fig. 1 B). 
Function-based classification of the identified proteins resulted 
in three major functional categories: (1) protein processing, 
(2) energy metabolism, and (3) oxidative stress, with some overlap 
between the latter two (Fig. 1 C). 

Biochemical evidence of strain-specific and 
□ DC-related changes in the expression of 
regulators of protein processing 

The proteomic data were subsequently subjected to further 
biochemical validation (Fig. 2). The expression of proteasome 
26S ATPase subunit 6 (PSMC6), which catalyzes the ATP- 
dependent 26S proteasomal degradation of ubiquitinated pro- 
teins (Smith et al., 2006), was significantly different between 
control and DDC-treated livers with respect to the existence of 
oligomeric PSMC6 complexes (Fig. 2 A, brackets), which is 
consistent with proteasomal impairments in this model (Harada 
et al., 2008), although there were no significant strain differences. 
At the mRNA level, PSMC6 was higher in the C3H strain under 
control conditions and elevated upon DDC treatment, whereas 
it was diminished upon DDC treatment in the C57BL livers 
(Table SI). The expression levels of monomeric protein disul- 
fide isomerase A4 (PDIA4), which catalyzes oxidative protein 



folding in the ER (Gruber et al., 2006), were identical between 
the two strains. In C57BL livers, there were slightly higher 
levels of selenium-binding protein 1 (SBP1; Fig. 2 A), which 
has been implicated in intra-Golgi transport and in selenium- 
dependent ubiquitination- and deubiquitination-mediated pro- 
tein degradation pathways (Porat et al., 2000; Jeong et al., 
2009). This was in contrast to the mRNA levels, which showed 
the opposite patterns (i.e., higher SBP1 mRNA in C3H livers) 
both before and after DDC treatment. The discordance between 
the protein and mRNA levels may reflect extensive posttransla- 
tional regulation of SBP1 and/or antibody reactivity with other 
SBP isoforms such as SBP2 (Lanfear et al., 1993). Although 
there were no major strain differences with regards to the total 
expression levels of PSMC6, PDIA4, and SBP1, the require- 
ment of energy and reducing equivalents for their proper func- 
tion suggests that differences may exist at the functional level. 
Therefore, we next examined the expression of energy metabo- 
lism regulators. 

Biochemical evidence of strain-specific and 
□ DC-related changes in the expression of 
energy metabolism enzymes 

Major differences in the expression of NDPK-A and -B iso- 
forms were observed between the two strains, with a significantly 
lower expression in C57BL livers under both control and DDC- 
treatment conditions (Fig. 2 B). The NDPK protein differences 
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Figure 3. Hepatic differences in the expression, activity, and distribution of the cytoplasmic enzyme NDPK. (A) Relative expression level of NDPK-B in 
extrahepatic organs from C3H and C57BL mice. The actin blot serves as a loading control. (B) Total NDPK activity in C3H and C57BL mouse livers under 
basal conditions and after DDC treatment (top) and in untreated isolated hepatocytes (bottom]. **, P < 0.01 using an unpaired f test; ***, P < 0.001 using 
a two-way analysis of variance. Each tested group included three to four mice, and samples were analyzed in triplicates. Results are represented as the 
mean and the SD. (C) Primary cultured hepatocytes from C3H (panels a-c) and C57BL (panels d-f) mice were triple immunostained for NDPK-B and K8 
and were then mounted in the presence of DAPI as described in Materials and methods. Bar, 10 urn. 



were not reflected at the mRNA level (Table SI), implying a 
posttranslational control of NDPK expression. NDPK activity 
is important for the synthesis of non-ATP nucleoside triphos- 
phates (i.e., CTP, GTP, and UTP) via the transfer of a phos- 
phate group from ATP (Boissan et al., 2009). NDPK also plays 
a role in the protection from oxidative stress-induced damage 
(Arnaud-Dabernat et al, 2004; Lee et al., 2009). The striking 
strain differences in liver NDPK expression (Fig. 2 B) were liver 
specific, as kidney, lung, heart, and spleen NDPK-B levels were 
similar between the two strains (Fig. 3 A). Further, total NDPK 
activity was significantly higher in C3H livers after DDC treat- 
ment in comparison with control C3H and DDC-treated C57BL 
livers (Fig. 3 B, top). The differences were more dramatic when 
assessing total NDPK activity in isolated C57BL hepatocytes, 
which was threefold lower compared with C3H hepatocytes 
(Fig. 3 B, bottom). NDPK exhibited a cytoplasmic filamentous 
organization and partial colocalization with keratin filaments 
in the hepatocytes (Fig. 3 C), although coimmunoprecipitation 
experiments did not show a direct interaction between the two 
(not depicted). 

In addition to NDPK, the protein expression of the glyco- 
lytic enzyme GAPDH was higher in the control C3H mice than all 
other groups (Fig. 2 B). At the mRNA level, GAPDH was approxi- 
mately twofold higher at basal conditions in C3H livers versus 
C57BL and exhibited an induction in C3H and suppression in 



C57BL livers after DDC treatment (Table SI). In addition to its 
traditional role in glycolysis and ATP generation, GAPDH per- 
forms numerous other functions, including redox sensing, in 
which oxidative stress-induced aggregation and nuclear translo- 
cation of GAPDH mediate cell death (Nakajima et al., 2009). 
Given the multiple roles of the housekeeping enzymes NDPK 
and GAPDH across energy metabolism and oxidative stress re- 
sponses, these data highlight a potential cross-regulation between 
these two pathways in DDC-induced liver injury. 

Strain-specific and DDC-related changes in 
the expression of oxidative stress-related 
enzymes and reactive oxygen species CROSJ 

In both mouse strains, DDC significantly induced the expres- 
sion of the antioxidant GSTs Mu and Pi (GSTM and GSTP) 
enzymes, which are critical for the detoxification of xenobiotics 
(Fig. 2 C; Whalen and Boyer, 1998). These data, in combination 
with previous results showing mouse sex differences in the han- 
dling of oxidative stress (Hanada et al., 2010), support a central 
role for oxidative stress in DDC-induced liver injury. Other anti- 
oxidant enzymes, including carbonic anhydrase 3 (CA3) and per- 
oxiredoxin 6 (PRDX6), which are induced by oxidative stress and 
protect against hydrogen peroxide-induced damage (Raisanen 
et al, 1999; Wood et al., 2003; Chowdhury et al., 2009), had lower 
mRNA (Table SI) and protein (Fig. 2 C) expression in C57BL 
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Figure 4. Strain-specific differences in oxidative stress-related enzymes and ROS. (A, left) Equal amounts of total liver protein from three untreated and 
four DDC-fed C3H and C57BL mice were resolved by SDS-PAGE and immunoblotted with antibodies to Cyp2el or Cyp3a. The PRDX6 blot serves as a 
loading control, (right) Relative band intensity values of the blots were plotted for each strain/treatment group. ***, P < 0.001 using a two-way analysis 
of variance. Each tested group included three to four mice. Results are represented as the mean and the SD. (B) Isolated primary hepatocytes from C3H 
(panels a-d) and C57BL (panels e-h) mice were cultured in the presence of DDC or vehicle (0.1% DMSO), loaded with the ROS indicator CM-H 2 DCFDA 
(panels a, b, e, and f) or immunostained for CBR3 (panels c, d, g, and h), and then mounted in the presence of DAPI as described in Materials and methods. 
Qualitative assessment of ROS levels (arrows in b and f) is based on CM-H 2 DCFDA fluorescence (green). CBR3 appears as perinuclear particulates (arrows 
in panel h). Bars, 20 urn. 



livers. Furthermore, C3H and C57BL livers differed significantly 
with respect to the presence of charged isoforms of PRDX6 
(Fig. S 1 A), which is consistent with the known functional regu- 
lation of PRDX6 by phosphorylation (Wu et al., 2009) and car- 
bonylation (Wong et al., 2008). Of note, a major strain difference 
was observed in the expression of carbonyl reductase 3 (CBR3; 
Fig. 2 C), which catalyzes the reduction of carbonyls and is 
known to be regulated by the antioxidant nuclear factor eryth- 
roid 2-related factor 2 (Nrf2) signaling pathway (Ebert et al., 
2010; Klaassen and Reisman, 2010). A significant (eightfold) 
induction in liver CBR3 expression was noted only in the C57BL 
but not the C3H mouse livers. Consistent with the lack of CBR3 
induction, C3H livers exhibited increased protein carbonylation 
after DDC treatment, which was not the case in the C57BL 
livers (Fig. SI B). 

Next, we examined the expression of prooxidant liver en- 
zymes and ROS generation. The expression of Cyp2e 1 , a known 
ROS-generating enzyme implicated in the pathogenesis of 
ALD (Cederbaum et al., 2009), and Cyp3a, the enzyme respon- 
sible for mouse liver DDC metabolism (Hanada et al., 2010), 
was significantly different between the two strains after DDC 
treatment (Fig. 4 A). Specifically, Cyp2el expression was 
down-regulated by 90% in the C3H and by only 30% in the 



C57BL livers, resulting in a fivefold higher relative Cyp2el ex- 
pression in C57BL compared with C3H livers after DDC treat- 
ment. The apparent decrease in Cyp2el protein was not a result 
of the contribution of bile ductular epithelium to total liver 
protein, as induction of K19 (ductal marker) was similar be- 
tween the two strains (Fig. S2 A), which is in agreement with 
previous findings (Hanada et al., 2008). Cyp2el mRNA levels 
also decreased in both strains (Fig. S2 B). In C57BL but not 
C3H livers, there was a marked Cyp3a induction (Fig. 4 A), 
which we have previously shown to correlate with the propen- 
sity of male mice to form more MDBs compared with female 
mice (Hanada et al., 2010). The differences in ROS-generating 
cytochromes P450 (CYPs) and GST induction suggested possi- 
ble glutathione depletion in the C57BL livers. However, mea- 
surement of the total and oxidized glutathione levels between 
the two strains did not reveal major differences (Fig. S2 C). 

Next, we examined the intracellular ROS levels ex vivo in 
cultured hepatocytes using the probe CM-H 2 DCFDA (Emslanov 
and Kusmartsev, 2010). As shown in Fig. 4 B (panels a, b, e, 
and f), DDC triggered a dramatic and more robust increase in 
ROS levels in C57BL hepatocytes in comparison with C3H 
hepatocytes. The induction of ROS correlated proportionately 
with increased CBR3 immune staining in DDC-treated C57BL 
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Figure 5. Differences in GAPDH aggregation 
and Siahl levels in C57BL and C3H mouse livers 
and effect of pioglitazone on DDC-modulated 
GAPDH nuclear localization in C57BL hepa- 
tocytes. (A) Nuclear fractions were prepared 
from C3H and C57BL mouse livers (livers from 
control diet [n = 3] or DDC-fed mice [n = 4]; 
groups are separated by dotted lines) and ana- 
lyzed on the same gel by SDS-PAGE followed 
by immunoblotting for GAPDH under reducing or 
nonreducing conditions. Lamin Bl was used as a 
loading control. Significant levels of high molecu- 
lar weight nuclear GAPDH aggregates were de- 
tected only in C57BL livers after DDC exposure. 
(B) C57BL livers (three independent control livers/ 
strain and four separate livers from DDC-fed 
mice) express significantly higher levels of Siahl 
protein after DDC treatment. Coomassie stain 
serves as loading control. (C) Biochemical analy- 
sis on total, cytoplasmic, and nuclei-enriched frac- 
tions from C57BL hepatocytes. The hepatocytes 
were cultured in the presence of vehicle (DMSO), 
100 pM DDC, or 3 pM pioglitazone (Pio) plus 
DDC for 48 h. Detergent lysates were then pre- 
pared and blotted with antibodies to GAPDH, 
lamin Bl (nuclear marker), p-tubulin (cytoplasmic 
marker), and pan-actin (loading control). 
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hepatocytes (Fig. 4 B, panels c, g, d, and h), which is consistent 
with the immunoblotting data (Fig. 2 C) and the dramatic in- 
crease in CBR3 mRNA levels (Table SI). These findings indi- 
cate that there is elevated oxidative stress in C57BL hepatocytes, 
triggering activation of antioxidant responses. 

Presence of nuclear GAPDH aggregates 
in DDC-treated C57BL mouse livers and 
isolated hepatocytes 

Given the major oxidative and metabolic differences between 
the MDB-resistant and MDB-susceptible mice, we focused our 
attention on the status of GAPDH as a potential common regu- 
lator of these pathways. Prior studies showed that oxidative stress 
induces extensive GAPDH modification and aggregation, which 
leads to its inactivation and compromised glycolysis (Nakajima 
et al., 2007, 2009; Tristan et al., 201 1). Furthermore, upon nitro- 
sylation (Hara et al., 2005) and association with the E3 ubiquitin 
ligase Siahl, GAPDH translocates to the nucleus and regulates 
autophagy induction (Colell et al., 2007), nuclear protein nitro- 
sylation (Kornberg et al., 2010), and apoptotic cell death (Hara 
et al., 2005). To determine whether GAPDH aggregation and 
nuclear translocation occur in vivo after DDC treatment, we ana- 
lyzed the nuclear fractions of livers from both mouse strains 
and found that total and aggregated nuclear GAPDH were sig- 
nificantly more abundant in C57BL livers (Fig. 5 A). Of note, 
the GAPDH aggregates were observed only when analyzed 
under nonreducing conditions, which is consistent with their 
intermolecular disulfide-bonded nature. The increased C57BL 
nuclear GAPDH after DDC exposure (Fig. 5 A) correlated with 
significantly higher levels of Siahl protein (which is involved 
in GAPDH nuclear translocation) in C57BL livers after DDC 
treatment (Fig. 5 B). Biochemical analysis of C57BL hepatocyte 
total, cytoplasmic, and nuclei-enriched fractions showed that 
DDC caused a significant reduction in the total and cytoplasmic 
levels of GAPDH to 12 and 14%, respectively, as compared with 



vehicle levels (Fig. 5 C). In contrast, there was an increase of 
215% in the nuclear levels of GAPDH after exposure to DDC. 
Of note, cotreatment with the insulin sensitizer pioglitazone 
almost completely reversed the total and cytoplasmic DDC- 
induced GAPDH depletion to 80% of vehicle and decreased the 
nuclear accumulation to 23% of vehicle (Fig. 5 C). Collectively, 
these data demonstrate a pharmacologically amenable role for 
GAPDH in hepatocellular injury. 

We also investigated the presence of GAPDH in the 
detergent-insoluble fraction of isolated hepatocytes and 
compared it with actin in the same fractions. C57BL, but not 
C3H, hepatocytes contained significant amounts of insoluble 
GAPDH (Fig. 6 A). Furthermore, immunofluorescence stain- 
ing showed that DDC induced the nuclear translocation of 
GAPDH in C3H hepatocytes (Fig. 6 B, panels a, b, e, and f). 
In stark contrast, GAPDH is already present in the nuclei of 
C57BL hepatocytes under basal conditions (Fig. 6 B, panels 
c and d) but undergoes aggregation upon exposure to DDC 
(Fig. 6 B, panels g and h). 

Strain-dependent effect of DDC on ex vivo 
hepatocyte expression and activity of 
energy metabolism and antioxidant enzymes 

We performed additional ex vivo experiments of isolated he- 
patocytes from C3H and C57BL mice to further explore the 
effects that were seen after DDC treatment in vivo. Similar to the 
in vivo findings (Fig. 2), ex vivo analysis showed that NDPK-B, 
CA3, and PRDX6 expression decreased after DDC treatment, 
whereas PDIA4 levels remained unchanged (Fig. 7 A; quan- 
tification is shown in Fig. S3 A). There was also a significant 
dose-dependent decrease in GAPDH and the tyrosine catabo- 
lism enzyme fumarylacetoacetate hydrolase (FAH) expression 
upon DDC treatment (Fig. 7 A), indicating a significant impair- 
ment in metabolic activity. These changes were more striking 
in C57BL hepatocytes, particularly in the case of NDPK-B. 
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Figure 6. Presence of detergent-insoluble and nuclear GAPDH in C57BL hepatocytes. (A) Hepatocytes from C57BL and C3H mice were cultured in pres- 
ence of the indicated concentrations of DDC or 0.1% DMSO vehicle (0) or in the absence of DDC/DMSO (Control). Detergent-insoluble equal fractions 
were prepared and then blotted with antibodies to actin (loading control) and GAPDH. (B) Immunofluorescence-based localization of GAPDH in vehicle 
(DMSO; panels a-d) or DDC-treated hepatocytes (panels e-h) that were also double stained with DAPI (nuclear staining). Note that after DDC treatment, 
GAPDH translocates to the nucleus in C3H hepatocytes (arrows in panels e and f) and appears to aggregate in the perinuclear region of C57BL hepato- 
cytes (arrows in panels g and h). Bar, 20 (jm. 



DMSO vehicle alone strongly decreased GAPDH, NDPK-B, 
and CA3 levels in C57BL hepatocytes (Fig. 7 A). The mecha- 
nism behind the DMSO effect is unclear; however, DMSO is 
known to exert control over the expression of several nuclear 
transcription factors (Su and Waxman, 2004), which may have 
contributed to its effects observed herein. DDC treatment aug- 
mented the expression of Cyp2el in both C3H and C57BL 
hepatocytes (Fig. 7 A). 

GAPDH is an upstream regulator of 
metabolic and antioxidant enzymes and 
RDS formation in hepatocytes 

We tested whether GAPDH is an upstream regulator of the 
proteins detected in our proteomic analysis (Figs. 1 and 2), 
given the known regulatory functions for GAPDH and the 
major differences in insoluble and nuclear GAPDH observed 
herein. As shown in Fig. 7 B (compare the first and second 
lanes with the fourth and fifth lanes), siRNA-mediated knock- 
down of GAPDH in isolated C57BL hepatocytes caused a sig- 
nificant down-regulation in the expression of NDPK-B (as well 
as NDPK-A; not depicted), FAH, CA3, and PRDX6 but had no 
effect on the expression of PDIA4 or Cyp2el when compared 
with control siRNA-transfected hepatocytes (quantification 
is shown in Fig. S3 B). Notably, the knockdown of GAPDH 
mimicked the effects of DDC treatment, thereby implicating 



GAPDH as an upstream regulator of the metabolic and anti- 
oxidant DDC-induced responses in hepatocytes. Consistent 
with this, GAPDH knockdown and overexpression increased 
2.2-fold and decreased 10-fold, respectively, the generation of 
ROS in DDC-treated C57BL hepatocytes (Fig. 7, C-E). 

□ DC-induced ROS formation in C3H and 
C57BL hepatocytes is augmented by 
GAPDH and NDPK knockdown 

We compared ROS formation in DDC-treated C3H and C57BL 
hepatocytes upon NDPK and GAPDH knockdown. The re- 
maining protein levels were 1.5-5% for GAPDH and 23-55% 
for NDPK after the treatment with the specific siRNA (Fig. 8 A). 
NDPK knockdown also caused a decrease in GAPDH levels 
to 68% in C3H and 41% in C57BL hepatocytes, indicating 
a coregulation between GAPDH and NDPK, which may be 
a result of changes in the cellular availability of nucleoside 
triphosphates. For example, a decrease in NDPK, which uses 
ATP to generate non-ATP nucleoside triphosphates (Boissan 
et al., 2009), could result in an increase in cytoplasmic ATP, 
which is an allosteric inhibitor of rate-limiting enzymes in 
glycolysis and thus may affect GAPDH levels (Berg et al., 
2002). The representative images in Fig. 8 B show elevated 
ROS levels after GAPDH or NDPK knockdown in DDC- 
treated C3H and C57BL hepatocytes. Quantification of the 
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Figure 7. DDC-induced ROS formation is associated 
with a depletion of energy metabolism and antioxidant 
enzymes and is regulated by GAPDH. (A) Hepatocytes 
from C3H and C57BL mice were left untreated (— ) 
or treated with either 0.1% DMSO vehicle (0) or the 
indicated concentrations of DDC for 48 h. Equal pro- 
tein amounts (NP-40 lysates) were analyzed for ex- 
pression of the indicated proteins. Samples from each 
strain were analyzed on separate gels. PDIA4 levels 
were unaltered (loading control). (B) C57BL hepato- 
cytes were transfected with Control or GAPDH siRNA 
for 24 h and were then cultured for an additional 
24 h in the presence of 0.025% DMSO vehicle (0) 
or 100 uM DDC or were left untreated ( — ). The 
NP-40 cell lysates were analyzed for the expression of 
the various proteins indicated. (C] Overexpression of 
Flag-tagged mouse GAPDH in isolated hepatocytes. 
NP-40 lysates were analyzed for expression of Flag- 
tagged GAPDH and total GAPDH. The Coomassie 
stain is included as a loading control. (D) Flag-GAPDH 
expression in C57BL hepatocytes, as determined by 
immunostaining with a mouse anti-Flag antibody (rep- 
resenting overexpressed GAPDH; green) and a rab- 
bit anti-GAPDH antibody (representing total GAPDH; 
red). (E) C57BL hepatocytes were mock transfected 
(Control) or transfected with GAPDH siRNA or Flag- 
GAPDH mouse cDNA for 24 h and were then treated 
with 1 00 uM DDC for an additional 24 h. Representa- 
tive images of ROS signal (green) with DAPI nuclear 
counterstain (blue) are shown. ROS levels (quantified 
as described in Materials and methods) exhibited a 
2.2-fold increase after GAPDH knockdown and a 10- 
fold decrease after GAPDH overexpression relative to 
control. Bars, 20 urn. 
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data from three independent images of each treatment group 
revealed statistically significant increases in ROS levels in C3H 
and C57BL hepatocytes upon GAPDH and NDPK knock- 
down (Fig. 8 C). 

GAPDH aggregates and NDPK complexes 
are present: in cirrhotic but: not: normal 
human livers 

To examine the potential clinical relevance of our findings, we 
compared the presence of cytoplasmic and nuclear GAPDH 
aggregates in livers from alcoholic cirrhosis patients with nor- 
mal control livers. Notably, the diseased but not normal human 
livers contained cytoplasmic and nuclear GAPDH aggregates 
(Fig. 9 A). Although there were no significant differences in 
monomeric NDPK-B protein between normal and diseased livers, 



the latter contained significant levels of higher molecular weight 
NDPK-B-containing species (Fig. 9 B). Future identification 
of the components of these complexes may yield additional 
mechanistic insight into a potential role of this enzyme in liver 
disease. Collectively, our findings in the genetically susceptible 
mouse model of MDB formation may extrapolate to human 
ALD with respect to GAPDH and NDPK function. 

Discussion 

In the present study, we tested the hypothesis that genetic fac- 
tors critically modulate susceptibility to hepatocyte injury in a 
mouse MDB model. We began with a global approach using 2D 
DIGE comparison of the basal protein expression levels and 
DDC-induced changes in livers of C3H (MDB resistant) and 
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Figure 8. DDC-induced ROS formation in C3H and 
C57BL hepatocytes is augmented by GAPDH and 
NDPK knockdown. (A) Biochemical analysis of protein 
expression of GAPDH, NDPK-B, and actin (loading 
control) 24 h after transfection of the isolated primary 
hepatocytes with the indicated siRNA. (B) Primary 
hepatocytes were transfected with the indicated siRNA 
for 24 h followed by the addition of 100 pM DDC 
for an additional 24 h. Representative differential in- 
terference contrast fluorescence images of ROS levels 
(green) in C3H and C57BL hepatocytes with DAPI 
counterstain (blue) are shown. Bar, 20 urn. (C) Quan- 
tification of hepatocyte ROS levels upon GAPDH and 
NDPK knockdown and DDC treatment. **, P < 0.01 ; 
***, P < 0.001 using a one-way analysis of variance 
and relative to the respective control group. ROS levels 
in control-transfected DDC-treated C57BL hepatocytes 
were significantly higher compared with control- 
transfected DDC-treated C3H hepatocytes (P = 0.0009 
using an unpaired t test). Results are represented as 
the mean and the SD (n = 15 cells/group), and the 
data for each tested group are representative of three 
independent hepatocyte isolations. 



C57BL (MDB susceptible) mouse strains. Next, the 2D DIGE- 
validated findings were used as a starting point to evaluate the 
status of the most likely candidate proteins and signaling path- 
ways ex vivo in hepatocytes and in vivo in mouse and human 
livers. Our findings indicate that MDB susceptibility is likely 
related to basal deficiencies in energy metabolism pathways 
that are intensified by elevated oxidative stress, ultimately re- 
sulting in compromised energy-dependent protein folding/ 
degradation responses (Fig. 10). We propose a central upstream 
role for GAPDH and its oxidative stress-induced aggregation 
and nuclear translocation in liver injury-induced hepatocyte in- 
clusion formation. The changes in GAPDH are unlikely to be 
unique to alcohol-related liver injury because analysis of liver 
explants from patients with other etiologies of end-stage liver 
disease also had accumulation of aggregated GAPDH (unpub- 
lished data). 

Imbalance between oxidative stress 
load and antioxidant defenses in 
MDB-susceptible livers 

We previously described an important role for microsomal 
CYP-mediated DDC metabolism and oxidative stress in the 
gender-dimorphic formation of mouse MDBs (Hanada et al., 
2010). Specifically, there is an induction of Cyp3a and sus- 
tained high expression of Cyp2el in livers of male mice, as 
contrasted to down-regulation of both of these enzymes in 
female mice upon chronic DDC exposure (Hanada et al., 2010). 



This results in increased Cyp3A-mediated DDC metabolism, 
oxidative stress (e.g., elevated lipid hydroperoxides), and a 
greater number of MDBs in the male mice (Hanada et al., 
2010). The current study lends further support to the role of the 
microsomal CYPs in the DDC-induced liver injury, as it relates 
to genetic background. 

Contribution by the other proteins that were identi- 
fied by the 2D DIGE analysis may have also significantly 
affected the histopathological changes that occur in the C3H and 
C57BL mice after DDC treatment. For example, exposure of 
C57BL hepatocytes to DDC in vitro resulted in a more dramatic 
decrease in the expression of several antioxidant enzymes, in- 
cluding PRDX6 and CA3, when compared with C3H hepato- 
cytes. This suggests that C57BL hepatocytes may experience 
higher H 2 0 2 levels as a result of deficient PRDX6 (Wood et al., 
2003) and H 2 0 2 -induced damage as a result of deficient CA3 
(Raisanen et al., 1999). Furthermore, expression of CBR3, 
a target of the Nrf2 antioxidant system (Ebert et al., 2010; 
Klaassen and Reisman, 2010), increased dramatically in DDC- 
treated hepatocytes and livers of C57BL mice in response to 
the elevated oxidative stress. Consistent with the CBR3 ex- 
pression differences, total liver protein carbonylation after 
DDC increased only in the DDC-treated C3H mice. Protein 
carbonylation plays important roles in cell signaling (Wong 
et al., 2010). For example, carbonylation of annexin Al targets 
it for proteasomal degradation, which in turn promotes cell 
survival and growth (Wong et al., 2008). 
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Figure 9. GAPDH and NDPK-B aggregates are present in nuclei-enriched 
fractions of human alcoholic cirrhosis livers. (A) Cytoplasmic and nuclei- 
enriched fractions from two different normal human livers (NL) and liver 
explants of eight patients with alcoholic cirrhosis were analyzed by SDS- 
PAGE followed by blotting for GAPDH (reducing [R] or nonreducing [NR] 
conditions). Lamin Bl and (3-tubulin were used as loading controls for the 
nuclear and cytoplasmic fractions, respectively. GAPDH aggregates (de- 
noted by asterisks) are present in the diseased but not the normal human 
livers. (B) Analysis of NDPK-B expression in total liver lysates (under reduc- 
ing conditions) from the same samples as in A. Higher molecular weight 
NDPK-B-containing complexes are present in the diseased but not the nor- 
mal human livers. 



response to oxidative stress, which is attributed to their ability to 
induce c-myc expression (Arnaud-Dabernat et al., 2004). We dem- 
onstrate herein that C57BL livers and hepatocytes exhibit a basal 
deficiency in NDPK, which is further accelerated by the presence 
of DDC, and that this contributes to elevated oxidative stress. Fur- 
thermore, multiple NDPK-containing protein complexes are pres- 
ent in cirrhotic but not normal human livers, implicating NDPK as 
a potential target in liver disease. 

Stress-induced GAPDH aggregation and 
nuclear translocation in MOB -susceptible 
mouse livers 

Aside from its role in glycolysis, GAPDH has numerous other 
important cellular functions, including regulation of cell sur- 
vival and cell death by way of its aggregation and nuclear 
translocation (Sen et al., 2008, 2009). Under conditions of 
oxidative stress, including ethanol-induced stress (Ou et al., 
2010), GAPDH translocates to the nucleus or forms amyloid- 
like aggregates in an irreversible process involving an active 
site cysteine residue (Nakajima et al., 2009). Neuroprotective 
agents, such as selegiline and rasagiline, block the ethanol- 
induced nuclear translocation of GAPDH, thus preventing the cell 
death-promoting signaling cascade (Ou et al., 2010). This cas- 
cade involves the formation of a complex between S-nitrosylated 
GAPDH with an NLS -containing E3 ubiquitin ligase (Siahl), 
translocation of the complex to the nucleus, and activation of 
p300/CREB-binding protein and its target, p53, to induce cell 
death (Sen et al., 2008). 

Our studies in the DDC-injured livers and isolated hepa- 
tocytes demonstrate that the redox-sensing role of GAPDH, 
which has thus far been described in neuronal systems, extends 
to hepatocytes in the liver. We find that depletion and nuclear 
accumulation of GAPDH are reversed by pioglitazone, which 
may provide a potential mechanistic basis for the histological 
improvements in livers of NASH patients who have been treated 
with pioglitazone (Promrat et al., 2004; Aithal et al., 2008; Sanyal 
et al., 2010). For example, in the PIVENS (pioglitazone versus 
vitamin E versus placebo for the treatment of nondiabetic 
patients with NASH) trial, a large placebo-controlled clinical 
trial of pioglitazone in NASH (Sanyal et al., 2010), there were 
significant improvements in steatosis, lobular inflammation, 
and hepatocellular ballooning in the pioglitazone compared 
with the placebo group (P < 0.001), which represents a potential 
clinical context for our findings. 



MOB susceptibility parallels liver- and 
hepatocyte-specific differences in NDPK 

NDPKs are structurally and functionally conserved enzymes from 
bacteria to humans, and they regulate intracellular energy pools by 
catalyzing the transfer of ■v-phosphates of a nucleoside triphosphate 
donor onto a nucleoside diphosphate acceptor and generate CTP, 
GTP, and UTP (Boissan et al., 2009). Aside from this homeostatic 
function, NDPKs are known to regulate transcription (Postel et al., 
1993) and cell proliferation (Braun et al., 2007) and are a novel 
therapeutic target for cancer metastasis (Steeg and Theodorescu, 
2008). Additionally, NDPKs exert protective effects in the cellular 



Materials and methods 

Antibodies 

The antibodies used in the study were to NDPK-A (NME1), NDPK-B 
(NME2), PSMC6, GAPDH (6C5), CYP3A1, FAH, PDIA4 (ERp72), lamin Bl, 
PRDX6, and Siahl (Abeam); GSTM1, CBR3, CA3, and ME1 (Santa Cruz 
Biotechnology, Inc.); panactin and (3-tubulin (Thermo Fisher Scientific); 
SPB1 (MBL International); Troma I (K8) and Troma III (K19; Developmental 
Studies Hybridoma Bank); mouse anti-DDK (OriGene); and CYP2E1 (a gift 
from H. Gelboin, National Institutes of Health, Bethesda, MD). 

Human and animal liver experiments 

The diseased human liver samples were explant tissues from patients 
who underwent liver transplantation for end-stage liver disease as a 
result of alcoholic cirrhosis and were used under an approved Human 
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Figure 1 0. Compromised energetics and redox cycling con- 
tribute to MDB formation in C57BL mice. DDC causes oxidative 
stress in vivo, as indicated by the significant induction of liver 
GST and CBR3 expression, and ex vivo in isolated hepato- 
cytes, as indicated by the accumulation of intracellular ROS. 
Elevated ROS promotes aggregation and nuclear translocation 
of GAPDH in the livers of DDGfed mice and in DDC-treated 
isolated hepatocytes. Ex vivo DDC treatment of hepatocytes 
causes a significant down-regulation of GAPDH. Furthermore, 
GAPDH is a central regulator of energy metabolism and oxida- 
tive stress-related enzymes, including NDPK, FAH, CA3, and 
PRDX6, as GAPDH knockdown leads to significant decreases in 
the expression of these proteins and mimics the effect of DDC. 
Therefore, diminished GAPDH expression and function lead to a 
bioenergetic and redox crisis, compromising energy-dependent 
protein processing pathways and ultimately leading to MDB for- 
mation. The glycolysis-stimulating agent pioglitazone reverses 
the GAPDH down-regulation and nuclear translocation. 



Subjects protocol (Ku et al., 2005). Nondiseased (normal) human liver 
tissues were obtained from the National Disease Research Interchange. 
Animal use was approved by and performed in accordance with the 
University Committee on Use and Care of Animals at the University of 
Michigan. DDC feeding for 3 mo to induce MDBs was performed as pre- 
viously described (Hanada et al., 2008). In brief, eight mice per strain 
were fed a powdered chow diet (Formulab Diet 5008; Dean's Animal 
Feeds) supplemented with 0.1% DDC wt/wt (Sigma-Aldrich) for 3 mo. 
Age-matched males (three mice per strain) were kept on a standard 
mouse diet (Teklad Global Diet 2019; Harlan Laboratories, Inc.) and 
were used as controls. 

2D DIGE and MS-based proteomic comparison of C3H and C57BL livers 

2D DIGE was performed as previously described (Marouga et al., 2005). 
Liver homogenates from untreated and DDC-treated age-matched C3H and 
C57BL male mice were labeled with the fluorescent dyes Cy3 (green) and 
Cy5 (red), respectively. The samples were mixed and resolved by 2D gel 
electrophoresis, and the differentially expressed proteins (indicated by 
green or red fluorescence) were selected for identification by MS/MS 
analysis. Approximately 140 spots (circled in Fig. 1 A) were selected for 
MS-based identification. 

Hepatocyte isolation and treatments 

Primary hepatocytes from 12-14-wk-old male C3H and C57BL mice 
were isolated as previously described (Snider et al., 201 1 ). In brief, the 
mice were anesthetized with 50 mg/kg Nembutal (Lundbeck Inc.), and 
the liver was perfused (7 ml/min) with 30 ml of solution I (Hank's bal- 
anced salt solution containing 0.5 mM EGTA, 5.5 mM glucose, and 1 % 
penicillin-streptomycin) followed by 25 ml of solution II (Hank's bal- 
anced salt solution containing 1.5 mM CaCI 2 , 5.5 mM glucose, 1% 
penicillin-streptomycin, and 2,000 U of collagenase IV [Worthington 
Biochemical Corporation]). Upon gentle mechanical dispersion of the 
liver, the cell suspension was filtered through a 70-pm cell strainer and 
pelleted (500 rpm for 2 min at 4°C), and the cells were cultured in 
Williams' medium E containing 10% FBS and 1% penicillin-streptomycin 
on collagen-coated dishes. The hepatocytes were allowed to attach for 
8 h (at 37°C in 5% C0 2 ) before the addition of siRNA constructs, DDC 
(in DMSO), DMSO, or DDC plus pioglitazone. For GAPDH and NDPK 
knockdown, hepatocytes were grown in 6-well plates (biochemical ex- 
periments) or 4-well chamber slides (immunofluorescence experiments) 
and were transfected with either control or GAPDH or NDPK-A + NDPK- 
B siRNA (Santa Cruz Biotechnology, Inc.) using Lipofectamine RNAiMAX 
reagent (Invitrogen) with 30 pmol siRNA and 5 pi Lipofectamine (per 
well on the 6-well plate). The medium was replaced after 6 h, and the 
cells were incubated for an additional 1 8 h before the addition of DDC. 
Transfection of Flag-tagged mouse GAPDH cDNA (OriGene) was per- 
formed using Lipofectamine LTX (Invitrogen) according to the manufac- 
turer's instructions. siRNA and DNA-transfected hepatocytes were 
exposed to 100 uM DDC for 24 h. 

Preparation of liver and hepatocyte lysates and immunoblotting 

Livers were homogenized in ice-cold NP-40 buffer (150 mM sodium chlo- 
ride, 1% NP-40, and 50 mM Tris, pH 8.0) supplemented with protease 
inhibitors. The NP-40 soluble and insoluble fractions were separated by 
centrifugation at 14,000 rpm for 20 min (at 4°C). Isolated hepatocytes 



were either lysed in NP-40 buffer or subjected to subcellular fractionation 
using the NE-PER cytoplasmic/nuclear fractionation kit (Thermo Fisher Sci- 
entific) to obtain total, cytoplasmic, and nuclear-enriched fractions. Liver or 
hepatocyte lysates were resolved on 4-20% gradient SDS-PAGE gels and 
were then transferred onto polyvinylidene difluoride membranes, which were 
subsequently blocked and incubated with the designated antibodies. 

Extraction of RNA and quantitative real-time PCR 

RNeasy kits (QIAGEN) were used to extract the total RNA from the mouse 
livers. RNA was translated into cDNA using the TaqMan reverse transcrip- 
tion kit (Applied Biosystems). cDNA was then subjected to quantitative 
real-time PCR (MyiQ real-time PCR detection system; Bio-Rad Laboratories) 
and was amplified by Brilliant SYBR green master mix using gene-specific 
primers (Table S2). 1 8S RNA was amplified as an internal control. Be- 
cause the amplification efficiencies for the gene of interest and the internal 
control were approximately equal, the quantification was expressed as a 
ratio of 2A a/9ene to 2A c,/inlemd conlro1 , in which ACt/gene and ACt/internal 
control represent the difference between the threshold cycle of amplifica- 
tion for the gene of interest and internal control, respectively. Each tested 
group included three to four mice, and samples were analyzed in tripli- 
cates. Results are represented as the mean and the SD. 

Immunofluorescence staining and confocal imaging of primary hepatocytes 

After methanol fixation (for 10 min at — 20°C), the hepatocytes were air 
dried, and nonspecific binding was blocked by incubation in blocking 
buffer (PBS with 2.5% wt/vol BSA and 2% goat serum). Primary antibodies 
to GAPDH and CBR3 in blocking buffer were incubated with the cells for 
1 h (at 22°C) followed by three 5-min washes in PBS and then incubation 
with secondary Alexa Fluor-conjugated antibodies (for 30 min). After rins- 
ing in PBS, cells were mounted with ProLong Gold containing DAPI (Invit- 
rogen), imaged on a laser-scanning confocal microscope (FluoView 500; 
Olympus) with a 60x oil immersion (1.4 NA) objective, and magnified 
with FluoView software (version 5.0; Olympus). DAPI, Alexa Fluor 488, 
and Alexa Fluor 594 were excited with a 405-nm laser diode, 488-nm 
argon laser, and 543-nm HeNe green laser, respectively, and sequential 
scans were used to maximize signal separation. 

Glutathione, NDPK activity, ROS, and protein 
carbonylation measurements 

Total and oxidized glutathione levels were measured using a commercial 
kit (Cayman Chemical) via an enzymatic recycling method using gluta- 
thione reductase. NDPK activity in liver homogenates and hepatocyte lysates 
was measured by continuous spectrophotometric rate determination using 
a three-step reaction with ATP, thymidine 5'-diphosphate, phosphoenol- 
pyruvate, pyruvate kinase, and p-NADH. The extinction coefficient of 
(3-NADH at 340 nm (6.22) was used to calculate kinase activity. For ROS 
measurement, untreated and DDC-treated hepatocytes were loaded with 
10 uM CM-H 2 DCFDA (Invitrogen) for 20 min, washed, and then mounted in 
ProLong Gold with DAPI and immediately imaged by confocal microscopy 
as described in the previous section. Protein carbonylation was measured 
using a commercial kit (Cell Biolabs, Inc.). 

Data analysis 

The graph data were presented and statistically analyzed using Prism 5 
software (GraphPad Software). Photoshop (CS2; Adobe) was used to 
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quantify ROS by measuring the green channel (ROS) signal intensity 
from 15 cells of each treatment condition and dividing by the image 
areas (in pixels]. 

Online supplemental material 

Fig. SI shows strain differences in PRDX6 modification and total protein 
carbonylation in C3H and C57BL livers. Fig. S2 shows a strain comparison 
of Cyp2el levels, ductal proliferation, and glutathione levels. Fig. S3 
shows quantification of the immunoblotting data shown in Fig. 7. Table SI 
shows relative mRNA expression of selected genes in untreated and DDC- 
treated C3H and C57BL mice. Table S2 lists the primers used for quantita- 
tive real-time PCR. Online supplemental material is available at http://www 
.icb.org/cgi/content/full/jcb.201 1021 42/DC1 . 
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